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' Abstract 

5^ ' One method of studying the asymptotic structure of spacetime is to apply Penrose's confor- 

^ l' mal rescaling technique. In this setting, the rescaled Einstein equations for the metric and the 

, conformal factor in the unphysical spacetime degenerate where the conformal factor vanishes, 

^ . namely at the boundary representing null infinity. This problem can be avoided by means of 

' a technique of H. Friedrich, which replaces the Einstein equations in the unphysical spacetime 

by an equivalent system of equations which is regular at the boundary. The initial value prob- 
lem for these equations produces a system of constraint equations known as the conformal 
^SJ ■ constraint equations. This work describes some of the properties of the conformal constraint 

^ ■ equations and develops a perturbative method of generating solutions near Euclidean space 

' under certain simplifying assumptions. 

o 

1 Introduction 

A model for the asymptotic structure of spacetime was suggested by Roger Penrose in p3| (see 



X 



also 1 15 for a review of the development of these ideas) using the technique of conformal rescaling. 



Since the reader is by now familiar with the details of the conformal rescaling construction, only 
I enough will be said here to fix the notation to be used in the remainder of this article. The 

bX)' object under study will consist of an asymptotically simple spacetime: that is, a physical spacetime 

consisting of a smooth, time- and space-orientable Lorentz manifold M with metric g satisfying 
the conditions: 



5—1 

' 1. M is diffeomorphic to the interior of an unphysical spacetime M, which is a smooth Lorentz 

manifold M with metric g that has a boundary dM; 

2. there is a smooth function : Af — > R, the conformal factor, which satisfies g = i^~^g on M 
(the pull-back by the diffeomorphism M int{M) has been suppressed for convenience); 

3. r2 = but dr2 7^ on dM; 

4. every null geodesic on M acquires a future and past endpoint on dM. 

Furthermore, {M,g) will be assumed to satisfy Einstein's vacuum equation Ric(g) =^ 0. 

The conformal boundary dM describes null infinity by virtue of condition (4), and asymptotic 
properties of the physical spacetime in null directions can be examined by studying the properties 
of the unphysical spacetime near its boundary. To this end, one could use the fact that, due to 
the conformal equivalence with the physical spacetime, the quantities Q and g must satisfy the 
conformally rescaled version of Einstein's equation, namely that Ric{Q,~^g) = 0. However, this 
equation has the drawback that it is degenerate near the boundary of M because there — > 0, and 
is thus not ideally suited for analytic investigations of the nature of the spacetime at null infinity. 
One possible means of avoiding this difficulty is to use a technique developed by Friedrich [pj| , 
which aims to describe the geometry of the unphysical spacetime by means of a new, yet fully 
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equivalent system of equations derived from the equation Ric(fl^^ g) = that is formaUy regular 
at the boundary of the unphysical spacetime. These equations involve g, Q and several additional 
quantities and are known as the conformal Einstein equations. 

As with Einstein's equations in the physical spacetime, it is possible to attempt to solve the 
conformal Einstein equations in the unphysical spacetime by means of an initial value formulation, 
where appropriate initial data are defined on a spacelike hypersurface Z in M and then evolved 
in time. Again as in the physical spacetime, the conformal equations induce certain constraint 
equations on the initial data; these equations are known as the conformal constraint equations and 
consist of a complicated system of coupled nonlinear differential equations for the induced metric h 
and second fundamental form X of Z, the conformal factor restricted to Z, and several additional 
quantities. A particular case of interest is when Z is asymptotically hyperboloidal, i. e. Z intersects 
dM transversely. In this case, the evolution of the boundary of Z forward in time produces the 
conformal boundary of the unphysical spacetime, and global questions concerning the existence of 
classes of spacetimes satisfying the definition of asymptotic simplicity can be addressed. See [^s) 
or |l^ for a review of these ideas. 

The purpose of this article is twofold. First, it is to introduce the conformal constraint equations 
and to investigate some of their properties, which will be done in Section 2. It will be found that, in 
a certain sense, they describe in a coupled way two mathematical problems — namely, the elliptic 
boundary value problem for the conformal factor O and the constraint problem arising from the 
Gauss-Codazzi equations of Z. Furthermore, a simple geometric assumption will be shown to lead 
to a special case of the equations in which the first problem does not appear and the second is in 
the forefront. In this special case, the full system of conformal constraint equations reduces to a 
much simpler and smaller system of equations that will be called the extended constraint equations 
because they will turn out to be equivalent to the usual vacuum Einstein constraint equations 
satisfied by the metric and second fundamental form of Z. (Tackling the boundary value problem 
is at present beyond the scope of this article but will be considered in the future.) 

The second goal of this article is to set up a perturbative approach for generating solutions of 
the extended constraint equations in the neighbourhood of a known solution, but only in the case 
of time-symmetric data — the more general case will be handled in another future paper . This 
task will be accomplished in Section 3 and the main theorem proved in this section appears on 



page 3.1, Because the extended constraint equations are equivalent to the usual constraint equa- 
tions, the Main Theorem can be interpreted as a new way of finding solutions of these equations, 
and furthermore, it will turn out to be one that is completely different from the 'classical' (i. e. 
Lichnerowicz-York) method of handling them. This issue will be discussed further in the Section 
3. 



2 The Conformal Constraint Equations 
2.1 Deriving the Equations 

Suppose {M,g,il) is an unphysical spacetime satisfying the assumptions of asymptotic simplicity 
and thus that the metric and conformal factor satisfy the rescaled version of Einstein's equation 

Ric{n-^g)=0. (1) 

This section sketches briefly how equation (|l|) for g and leads first to the conformal Einstein 
equations for g, Q and additional quantities, and then to the conformal constraint equations. Begin 
by expanding (|l|) to obtain 

on 2 „ „ ^ sv^nvx^ 
— IT^'"" ^ n ^M^"-^ ^ — 5^1^ , (2) 
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where are the components of the Ricci tensor in the unphysical spacetime, is the covariant 
derivative of the four-metric and □ is its D'Alembertian operator. Notice that, as it is written, 
equation (j^) contains terms with negative powers of which tend to infinity near the boundary 
of the unphysical spacetime. Alternatively, if the equation is multiplied through by 17^, then the 
principal parts of the differential operators acting on g and fJ, would tend to zero at the boundary. 
Either way, equation (H) degenerates near the boundary of the unphysical spacetime, making it an 
unwieldy choice for studying the geometry of the spacetime near null infinity. 

Helmut Friedrich has developed a procedure for obtaining a system of equations equivalent to 
the rescaled Einstein equations (|^) but that is formally regular at the boundary of the unphysical 
spacetime and thus avoids the problems outlined above. This work can be found in several papers, 
see for example . Friedrich's derivation proceeded in the following way. Let C^uXp be the Weyl 
tensor of the metric g and define the quantities 

Sfj.u\p — ^ ^C^i^xp (3) 
4 24 

The tensor S^uXp is smooth on dM because under the assumptions of asymptotic simplicity, Penrose 
has shown that C^^xp vanishes at the boundary of M (a further condition on the topology of 
dM — that it admit spherical sections — is also needed, and will be assumed to hold). Friedrich 
then found that the system of equations 

(4) 

VSpx.p = 
- VpQV^'fl = 

RpuXp — ^S^ijXp + gp[xLv\p — Lp_[x9iy]p J 

where Rp^Xp is the Riemann curvature tensor of the metric g, can be derived from (j^). This 
is done by rephrasing (H) in terms of the quantities (^), adjoining the Bianchi identity for the 
physical spacetime and the unphysical spacetime, and adjoining the well-known decomposition of 
the curvature tensor given by 



Rpi/Xp — C^iyXp + gpxLup ~ 9pp^vx + dvpLpx ^ 9vxL 



pp ■ 



Then the algebraic properties of the equations and the defined quantities allows them to be ma- 
nipulated into the regular ones listed in (||), which are known as the conformal Einstein equations. 
The equivalence of (|) to (||) is confirmed as follows. Suppose the quantities L, S and ij) as well 
as g and Vl satisfy (^). Then by algebra, it can be shown that the pair {g,^) satisfies (||) and that 
L, S and relate to f2 and the curvature quantities in the manner indicated in (|^). (The algebra 
is fairly straightforward: for instance, the last equation in (^) identifies L and S as components of 
the curvature tensor; then it is a matter of computation to recover equation (^) from the remaining 
five.) 

It is immediately clear that the equations in ^ are regular when = 0. Furthermore, not 
only do the conformal Einstein equations contain the rescaled vacuum Einstein equations, but they 
also contain the Bianchi identity for the curvature tensor, though expressed in the new unknowns. 
Thus one can consider to contain compatibility conditions since the Bianchi identity is in some 
sense a compatibility condition for the curvature tensor — meaning that the Bianchi identity 
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is a result of requiring second covariant derivatives to commute properly (this can best be seen 
explicitly by rewriting the curvature tensor in terms of the vector- valued connection 1-forms as in 
Q, whereby the Bianchi identity becomes an incarnation of the identity d'^ = satisfied by the 
exterior differential operator). 

Suppose now that Z is a spacelike hypersurface in M. The fact that the conformal Einstein 
equations constrain certain initial data on Z can be seen by performing a 3+1 splitting of spacetime 
near Z. Choose a frame Ea, a — 1,2, 3, for the tangent space of Z and complete this to a frame for 
the unphysical spacetime by adjoining the forward-pointing unit normal vector field n of Z. Use 
this frame to decompose the equations into components parallel and perpendicular to Z. The 
constraint equations induced by the conformal Einstein equations are those equations in which no 
second normal derivatives of g or f2, and no first normal derivatives of L, S or appear. The 
initial data are the unknown quantities which are found in these equations; they are: 

• the induced metric of Z, which will still be called g (no confusion will arise because the 
4-dimensional setting will not be considered further in the remainder of the this article), 

• the second fundamental form X of Z, 

• the function restricted to Z, 

• the normal derivative n(f2)|^, to be denoted E, 

• the tensors Lab = E'^E'^L^^ and La = n'^E'^L^^, 

• the tensors Sabc - n^E'^^E^EPS^^xp and Sab = nf^n^'E^EPSx^p., 

• and the function ip restricted to Z. 

The constraint equations arising from the 3-1-1 splitting are: 

VaVbfi = Y,Xab - ^Lab + ijjgab 
VaE = XlVc^l - nLa 
- -^'^Lba - ^La 
^aLbc ^ ^bLac = V'^57S'eca6 ~ ^Scab ~~ {^acLb ~ ^bcLa) 
VaLb — ^bLa ~ V^fJ 5eab + ^a-^bc — ^b^ac 

(5) 

Vac Q -wa Q ^ ' 

^abc — ^fj'-'ac ~ -^c^ab 

Vac T/ac o 

'->ab — '-'abc 

= 2nij + j:^ - ||vi7f 

Vc^ba — VfjXca — rtSabc + gabLc — gacLb 
Rab = ^Sab + Lab + -^L^gab — ^c-^ab + ^ca^b 

where V now denotes the covariant derivative operator on Z corresponding to its induced metric g 
and Rab is its Ricci curvature. These equations are known as the conformal constraint equations. 
The derivation of these equations will not be reproduced here — the reader is asked to consult 
|l6| for this material. However, it is fairly easy to recognize the origin of the various terms 
appearing above. For example, the first two equations arise as the tangential and tangential- 
normal components of the first equation of (^). Furthermore, and more importantly for the sequel, 
the last two equations arise as the Gauss and Codazzi equations applied to the decomposition of 
the curvature tensor given by the last equation of (Q). 
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The various tensor quantities that appear in (|^) possess certain symmetries as a result of their 
origin as components of the curvature tensor: Lab is symmetric; Sab is symmetric and trace-free; 
and Sabc is antisymmetric on its last two indices, satisfies the Jacobi symmetry Sabc+Scab+Sbca = 
and is trace-free on all its indices. (Tensors with these symmetries will appear often in the sequel. 
Tensors of rank three that are antisymmetric on their last two indices and satisfy the Jacobi 
symmetry will be called Jacobi tensors for short while those which are in addition trace-free will 
be called traceless Jacobi tensors.) Note that even though the tensor Sabcd — E^^Ej^ E^E^iSf^j^xp 
appears in the constraint equations, it is not a truly independent initial datum because, thanks to 
the symmetries of ^^^Ap, it can be written as Sabcd = ga[cSd]b - Sa[cgd]b- 

The system is clearly exceedingly complicated because it is quasi-linear and highly coupled. 
However, the advantage provided by is once again that it is formally regular at the boundary 
of Z. For the sake of comparison, recall the interior of Z can be viewed as a spacelike hypersurface 
of the physical spacetime, and as such, satisfies the usual Einstein constraint equations there. In 
other words, if its induced metric is denoted by g and its second fundamental form by X, then 



mab - ^bK = 

(6) 

where V is the covariant derivative operator of the metric g and R is its scalar curvature. These 
equations can be rephrased in terms of g, X and f2 in the unphysical spacetime by conformal 
transformation. The necessary transformation rules are that g = il^^g and also that X = fl^^X + 
T,n^^g (which can be found by conformally transforming the definition of the second fundamental 
form as the normal component of the covariant derivative restricted to Z). The resulting equations 
are 

n^{R + (x;j)2 - x'^^'Xab) + 4nAgn - 6|| vrjf + 4i7sx;; + 6t.^ 

n{\/aX'^ - VfcX;^) - 2VbS - 2X^ nablaail = , 

where S = n(ri)|^ and Ag is the Laplacian of the metric g. Once again, the principal parts of 
these equations contain factors of 51 and thus degenerate as — > near the boundary of Z. This 
behaviour does not arise in the conformal constraint equations. 

The conformal constraint equations listed in (^ are equivalent to the usual constraint equations 
(0) because if (g, X, ft, S) solves (0) and the subsidiary quantities S, S, L and ip are defined as 
indicated in (H) (e. g. the last equation defines ip; then the first equation defines the 2-tensor 
Lab, etc. ), then straightforward computation shows that the conformal constraint equations are 
satisfied; furthermore, if {g,X,il,J^, S, S, L,i/j) satisfies (||), then it can be shown that {g,X,fl,l^) 
satisfies (|^), and consequently, g and X, given by the transformation rules above, satisfy the 
usual constraint equations (^. These considerations thus suggest one method for constructing 
solutions of the conformal constraint equations: construct any solution (g, X) of the usual constraint 
equations using standard techniques, choose a conformal factor, perform the transformations to the 
unphysical spacetime and use the conformal constraint equations to define the subsidiary quantities 
in terms of (g, X). Then these new quantities satisfy the conformal constraint equations. 

Consequently, it is possible to assume the existence of initial data with well-defined asymptotic 
properties (essentially given by the transformation rules above) and study only the time evolution 
of the data according to the conformal Einstein equations (^). This is the idea behind the work 
of Friedrich in (extended in ||l9|), where the time evolution of suitably small initial data on 
an asymptotically hyperboloidal hypersurface was studied and a complete future development was 
found. The nature of the asymptotic structure of this class of solutions near null infinity, and in 
particular the relationship between the asymptotic structure of the solution and the asymptotic 
structure of the initial data, was then analyzed extensively by Andersson, Chrusciel and Friedrich 
in Ql (extended by Andersson and Chrusciel in |j), and was based on the rescaled Einstein 
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equations ^ and their constraints (|7|). However, the problem of the vanishing of the conformal 
factor near the boundary of the unphysical spacetime and the resultant degeneration of these equa- 
tions remains a part of the ACF methods. Thus they are not ideally suited for certain applications, 
in particular for implementing numerical studies of asymptotically hyperboloidal data near null 
infinity where the presence of negative powers of can cause computational codes to crash (see 
|]l5| for details). It is for this reason that new methods for solving (|^) directly, rather than through 
the usual constraint equations, must be developed. This question will begin to be tackled in the 
remainder of this article. 



2.2 Reduction to the Extended Constraint Equations 

The complexity of the conformal constraint equations makes it a daunting task to attempt to 
develop any methods for obtaining solutions of the equations in their full generality. However, a 
great deal of structure is contained within these equations, and the hope is that this structure can 
be exploited in the search for solutions. For instance, it is possible to disentangle in some sense 
the equations relating to the conformal factor and its associated boundary value problem from 
the equations related to the Gauss-Codazzi equations of Z by restricting to a special case of the 
equations. 

The special case that will be considered in the rest of this article is to assume that the conformal 
diffeomorphism between M and M is the identity, and consequently that the conformal factor is 
trivial (i. e. = 1) in the unphysical spacetime. This is somewhat of a strange simplification, 
because it requires that the spacetime M have empty boundary (since 1I~^(0) = dM)\ One 
would thus not find oneself in this special case in practice since the whole point of the conformal 
constraint equations is to study hyperboloidal initial data in a conformally rescaled spacetime 
that has a boundary at null infinity. Nevertheless, the simplification afforded by the assumption 
f2 = 1 is worthwhile to consider from a mathematical point of view because it will accomplish the 
disentanglement described above and allow the Gauss-Codazzi-type equations within the conformal 
constraint equations to be studied in isolation. 

To see this explicitly, one must substitute = 1 and S = (which is consistent with the 
assumption that 17 = 1 in spacetime since E = n(ri)|2 = where n is the forward-pointing unit 
normal of Z) into the equations (|^). One first sees that Lab, La and "0 are forced to vanish under 
this assumption, and then that the conformal constraint equations reduce to the following system 
of four coupled equations: 

-Rafc = Sab ~ ^c'^ab + ^a^cb 

VcXdf, — VfcXac = Sabc 

Van ya q -ya q 

iJabc — ^b'^ac ~ ^c'^ab 

Van \'ac n 

Jab — Jabc ■ 

Here, covariant derivatives are taken with respect to the induced metric gab of Z and Xab is the 
second fundamental form of Z. As before, the tensor Sab is symmetric and trace-free with respect 
to gab whereas the tensor Sabc is a traceless Jacobi tensor. These four quantities are the unknowns 
for which these equations must be solved. For reasons that will become apparent later on, it will 
be helpful work instead with the equivalent system obtained by replacing Sab and Sac in the third 
equation by Rab and Rac from the first equation. The system one obtains is actually just 

Rab = Sab ~ -^c-^afc + ^'a^cb 
VcXaf, — VbXjjc = Sabc 

(8) 

Vac -ya D -ya r? 

•Jabc — ^b^ac ~ -^c^ab 

Van vac c 

'Jab — ~^ Jabc : 
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because the terms cubic in X vanish. 

Notice that because of the symmetries of S and S, if the traces of the first two equations of 
are taken, then the usual constraint equations (^) resuh. Furthermore, if gab and Xab satisfy 
the usual constraint equations and one defines Sabc and Sat by the first two equations of (|^) 
respectively, then the remaining two equations follow by straightforward algebra and the Bianchi 
identity. Thus equations (^) are equivalent to the usual vacuum Einstein constraint equations and 
for this reason are called the extended constraint equations. 

2.3 Properties of the Extended Constraint Equations 

The extended constraint equations (j^) are clearly formally much simpler than the full system of 
conformal constraint equations. However, several essential features of the full equations remain. 
These features refer to the ellipticity properties of the various differential operators appearing in 
as well as to the compatibility conditions built into these equations. 

Ellipticity Properties 

One must consider the principal symbols of the operators that appear on the left hand sides of 
the extended constraint equations in order to understand their ellipticity properties. Begin with a 
definition of the symbol. Recall that if P : C°°(R", R^) ^ C°°(R",R*0 is a linear differential 
operator of order m with constant coefficients, then it can be expressed as 

1=1 ^ ' / 

where Pq is a differential operator of order less than or equal to m — 1 and the 5"^ are elements 
of R''^. The principal symbol of P is the family of linear maps given by 

ai-\ \-an—m / 

for any non-zero (^i, . . . , £,n) G R" and v G R^. Furthermore, the operator P is called underdeter- 
mined elliptic if the symbol is surjective for each non-zero ^, overdetermined elliptic if the symbol 
is injective for each non-zero ^ and simply elliptic if the symbol is bijective for each non-zero ^. 
An operator with non-constant coefRcients has a symbol at each point of the domain, while for a 
nonlinear operator, it is the linearization which has a symbol at each given u e C°°(R", R^). Such 
operators are overdetermined, underdetermined or elliptic if their symbols possess these properties 
uniformly. 

To understand the ellipticity properties of the conformal constraint equations, begin with the 
equation for the metric gab- It is quasi-linear in g, with highest-order terms given by 

^ cd f 9^9ad d^gbd 1 d^gab 1 d^gcd \ 
9ab^g yg^bQ^c-^ Q^ag^c 2 dx^dx"^ 2dx''dx^)' 

The linearization of this expression at a given metric is neither over- nor underdetermined elliptic, 
nor is it elliptic. However, it is well known that the Ricci curvature is degenerate as an operator 
on metrics because it is invariant under changes of coordinates of the metric, and that the Ricci 
curvature operator can be made formally elliptic by making an appropriate choice of coordinate 
gauge. The standard choice is to require that the metric be expressed in harmonic coordinates, 
which are defined by the requirement that the coordinate functions a;° are harmonic functions, 
i. e. that l^hx"' = for each a. (Since the metric itself depends on the coordinate functions, the 
requirement that the coordinates be harmonic is in fact a nonlinear condition. Nevertheless, the 
existence of such coordinates, defined outside sufficiently large balls in R'^ for any asymptotically 
flat metric, has been guaranteed by Bartnik in H.) 
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To show that the Ricci operator is ehiptic in harmonic coordinates, first note that a straightfor- 
ward calculation implies that the harmonic coordinate condition Agx" = for all a is equivalent 
to the condition g^'^T'^^ — for all a on the Christoffel symbols of g. Now set = g^'^^bc (a-nd 
also Fa = gas^"), and then recall that the components of the Ricci tensor satisfy 

^ah = -R^ + 2 (ra;6 + Tfj^a) (9) 

where are the components of the reduced Ricci operator defined by 

Rab^-l9'''9ab,rs+qiT). (10) 

In the expressions above, a comma denotes ordinary differentiation with respect to the coordinates, 
a semicolon denotes covariant differentiation (since is not a tensor, this is to be taken formally; 
i. e. Ta-t = Ta^b — TsF^j), and q(T) denotes a term that is quadratic in the components F". The 
reduced Ricci operator is clearly elliptic in g. Since F" = for all a in harmonic coordinates, 
Rab{g) — Rabid) thcse Coordinates, and thus the Ricci operator is elliptic in g when g satisfies 
the harmonic coordinate condition. 

The second equation in the extended constraint equations is linear in Xab and its left hand side 
defines a differential operator X^fc ^c^ab ~ ^b^ac from the space of symmetric tensors to the 
space of Jacobi tensors. (It can be easily verified that the left hand side of the first equation in ^ 
satisfies the relevant symmetries. However, it can also be verified that the left hand side is not a 
priori traceless on all its indices — this is only a requirement on the eventual solution since the 
left hand side is equated with a traceless Jacobi tensor.) The principal symbol of this operator is 

: ^ab I— > ic^ab — ^b^ac ■ 

By the following simple argument, one can show that has a one-dimensional kernel and is not 
surjective. 

Suppose first that a^{Xab) = for some non-zero ^. Since ^a'?'^ 7^ 0, one can write uniquely 
Xab = Xl^j^ + cS^a^b for somc c, where XjJ^ is trace-free. Substituting this expression for Xab yields 

6xL-^cX°, = o. (11) 

Taking the trace over a and b implies that C^X"^ = 0. Then, contracting with ^"^ gives C^Cc^afc = 0, 
or 0. Consequently, the kernel of the symbol is one-dimensional, and consists of tensors 

of the form c^aCfc- Next, since the space of symmetric 2-tensors is six-dimensional, the image 
of the symbol is five-dimensional. Now, the target space of Jacobi tensors is eight-dimensional 
because any Jacobi tensor can be decomposed as Tabc = s'^bc^ae + ^bgac — Acgab where Fae is a 
trace-free and symmetric tensor (accounting for five dimensions), Ab is a 1-form (accounting for 
the remaining three), and Sabc is the fully antisymmetric permutation symbol. The symbol can 
thus not be surjective. Note, however, that when it is restricted to trace-free tensors, the principal 
symbol is at least injective. Consequently, the first equation of (^) is overdetermined elliptic when 
restricted to the space of trace-free symmetric 2-tensors. 

The third and fourth equations in (^) are linear in Sabc and Sab respectively. It can be shown 
that the operators Sabc ^ '^"'Sabc and Sab V^S'af, are underdetermined elliptic by demonstrating 
that their principal symbols Sabc '—>' £.°'Sabc and Sab * ^""Sab are surjective maps from the space 
of symmetric, trace-free tensors onto the space of 1-forms and from the space of traceless Jacobi 
tensors onto the space of antisymmetric 2-tensors, respectively. These are fairly straightforward 
calculations and left to the reader. 

Compatibility Conditions 

As mentioned in Section 2.1, the conformal Einstein equations (^ contain the Bianchi identity, 
and was interpreted as being a compatibility condition for the other equations. Such compatibility 
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conditions are also to be found in the conformal constraint equations; the present goal is to exhibit 
this explicitly. Begin by considering the first and fourth equations in (H) . The Bianchi identity for 
the Ricci curvature is 

whereby the first equation of implies 

= V^{Sab - X^Xab + XlX,b) ^l^b{- (^c)' + X'^'Xac) 

= V^Sab - (X^S^ - X^) (V"X,„ - VaXl) - X-'^iVbXac " VcXab) 

= VSab — {X^S^ — X^)h"^ Suav + X°''^Sabc (12) 

using the second equation in ^ and its trace. By the symmetries of Babe, the middle term in ( p^ 
vanishes, leaving 

— V^Sab + X°''^Sabc , 

which is exactly the fourth equation of (@). 

The second and third equations of (|^) consist of a constraint equation and its compatibility 
condition as well, but in a different sense. Recall that what a compatibility condition should 
reflect is that second second covariant derivatives commute properly. Consider, then, the result of 
commuting the second covariant derivatives of the second equation of (^). Begin with 

Sabc — '^cXab ~ VfoX^c 

and compute 

e^'"^'VeSabc — 'ie'^'^'^'S/ cXab 

^e'''"'{VeycXab-VcyeXab) 
= £'^''^(-Reca''Xs6 + R^^^^^Xas) 

^ Reca^Xsb (13) 



since the symmetries of Rabcd imply that £°'^'^Rabcd — 0. Now substitute in (f3) the well-known 
decomposition of the curvature tensor in three dimensions, namely that 

Reca" = 9eaRl - K^ca + S^Rea " ffca-Re " ^^{geaSc ~ Kdca) , 

to obtain 

e^''^V,Sabc = 2£''\XlR,,. (14) 

Claim: equation ( p^ ) is exactly the third equation of (||). To sec this, recall that a traceless Jacobi 
tensor can be decomposed as Babe = £'^bc^o-e where Fae is trace-free and symmetric. Consequently, 

e^^^V eBabc = ^'^^'^'^ e^^bc^au 
= ^V'Fae 

= 2V^i^ea (by symmetry) 

c- fccyye JI Tp 

— V e bc^eu 

= Sa'^V^Bebc ■ (15) 

Thus (|l|) together with (|l|) implies that 

^a"^'^'^ Bcbc — '^s'"^ aXb^cs , 

which is the third equation of M) (or at least its dual, but this is equivalent). 
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3 Asymptotically Flat Solutions of the Extended Constraint 
Equations in the Time Symmetric Case 

3.1 Statement of the Main Theorem 

Because the conformal boundary of the spacetime M is absent under the triviahty assumptions 
that have been made on the conformal diffeomorphism, a natural setting in which to investigate the 
extended constraint equations (||) is the case in which M is asymptotically Minkowski space and 
that Z is asymptotically flat. In fact, one solution of the extended constraint equations satisfying 
these conditions is when Z = R"^ and the initial data is Euclidean metric g ~ S with vanishing 
second fundamental form and tensors S and S. Neighbouring asymptotically flat solutions are those 
whose metric 5 is a small perturbation of S that decays suitably to S near infinity, and X, S and S 
are also small and decay suitably. These solutions are in addition time symmetric if their second 
fundamental form actually vanishes identically. The theorem that will be proved in the remainder 
of this article is a characterization of the space of asymptotically flat and time-symmetric solutions 
of the extended constraint equations in the neighbourhood of the trivial solution given above. The 
case of non-time-symmetric solutions is as yet beyond the scope of this article, though a future 
paper by the Author will clear this up 0. 

Under the assumption of time-symmetry, the requirement that X = implies that 5 = as 
well, and so the extended constraint equations further reduce to the following system of equations: 

V'^Sab = 
Rab{g) — Sab 

for the unknown metric g and unknown trace- free and symmetric tensor S. Since these equations 
will be solved for metrics near the Euclidean metric, it will be preferable to write metrics as small 
perturbations of the Euclidean metric of the form 6 + h where /i is a symmetric tensor suitably 
near 0. Thus the above system should be replaced with the system 

V^'Sab = 

(16) 

Rab{S + h) ^ Sab ■ 

The covariant derivative here corresponds to the metric S + h. The theorem that will be proved is 
the following. 

Main Theorem: There exists a Banach space B of free data along with a neighbourhood U of 
zero in B, Banach spaces Y and Y' of symmetric 2-tensors, and smooth functions : U Y and 
ij)' : U Y' with "0(0) ~ ?A'(0) ~ Q so that for every b E U , the following hold: 

1. ^j{b) and ij/ib) tend asymptotically towards zero; 

2. g = S + 'ip{b) defines an asymptotically fiat Riemannian metric on 'Rp ; 

3. S = ip'{b) defines a symmetric tensor that is trace-free with respect to g; 
4- g and S satisfy the equations (p^). 

The proof of this theorem will be presented in the remaining sections of this article, and consists of 



essentially two steps. As outlined in the previous section, (16) is not an elliptic system. However, 
by exploiting the elliptic properties of the equations, it is possible to define a closely related system 
of equations, called the associated system, which is elliptic. In it, the tensor S is decomposed into 
a sum of two components of the form T + P{X), where T is a symmetric and trace-free tensor, X 
is a 1-form and P is the adjoint of the divergence operator Sab V^^ab. The system ([l^), written 
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in terms of this decomposition, yields equations for g, X, and T whose hnearization in the g and X 
directions is bijective (or near enough to being bijective — this wiU be cleared up in due course). 
Thus the Implicit Function Theorem can be invoked to find solutions where the quantities g and 
X are expressed as functions of T, which consists of the first step of the proof. The second step 
is then to show that all solutions of the associated system are also solutions of the original system 
([l6|). The Author wishes to thank H. Friedrich for suggesting this approach for solving ([l^). 

The method outlined above for solving the extended constraints in the time symmetric case is 
in fact a method for solving the usual vacuum constraint equations in the time-symmetric case 
(namely the equation R{g) = 0, which follows from (^ by taking a trace) because of the equiv- 
alence of the extended constraints and the usual constraints described earlier. The differences 
between this method and the 'classical' Lichnerowicz-York method for solving the constraint equa- 
tions are now readily apparent. In the classical method, one freely prescribes a metric go on R"^ and 
considers the conformally rescaled metric g = u^go, where u : Y{? ^ R is an unknown function. 
One then reads the equation R{u'^gQ) = as a semi-linear elliptic equation for u. In contrast, the 
present method treats the metric g and the one-form X as the unknowns and leads to a quasi- 
linear elliptic system for these quantities in terms of the freely prescribable quantity T, which is a 
component of the curvature of the solution. 

Remark: The Main Theorem does not fall into the domain of prescribed Ricci curvature as, for 
example, do the results of De Turck and his collaborators In these papers, the 

authors suppose a fixed symmetric tensor S is given on a set O and attempt to find conditions 
under which a metric g exists on O so that Ric{g) — S. In the Main Theorem, by contrast, the 
tensor S is itself an unknown quantity and only a component is prescribed ahead of time by the 
free data. Furthermore, De Turck's results are local in nature since O is usually an open set in 
R", while the Main Theorem gives a global (though perturbative) result. 



3.2 Formulating an Elliptic Problem 

The first task in the proof of the Main Theorem is to construct the associated elliptic system that 
is to be solved by the Implicit Function Theorem. What is needed is a system of equations closely 



related to (16) but that is elliptic. To this end, the the Ricci curvature operator in (|16| ) will be 
replaced by the reduced Ricci operator, which is elliptic as described in Section 2.3. Making this 
substitution is equivalent to assuming a priori that the harmonic coordinate condition is satisfied 
by the metric S + h. Of course, this assumption must be justified later on; i. e. it must be shown 
that 6 + h does indeed satisfy the harmonic coordinate condition, and this is the intent of the 
second step of the proof of the Main Theorem. The remaining operator in (^) is underdetcrmined 
elliptic, and an elliptic operator can be constructed from this by using a standard technique known 
as the York decomposition (see |^ but also ^ for a thorough analysis of this method) . Write 
a symmetric, trace-free tensor 5" in terms of a 1-form X and a freely prescribed tensor symmetric 
T as 

S{h,X, T) = T* +C^+''{X) . 

where T* = T - ^TTs+hiT){S + h) is the trace-free part of T and C^+''{X) is the conformal Killing 
operator with respect to the metric S + h acting on X. This is defined for a general metric g by 

where V is the covariant derivative of the metric g. The reason for making this choice is that 
the composition of the divergence operator in (jl^) and the conformal Killing operator, that is the 
composite operator divg o £f given componentwise by 

Xa ^ V°(VaXfc + VbXa - ^V^X,5„fc) = V'V.Xfc + ^VfeV^X, + Rl{g)X, , 
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is elliptic, as can easily be seen by computing its symbol or by making the following observation. 
It easy to compute that the conformal Killing operator is the formal adjoint of the divergence 
operator Sab '-^ V^S'ab taking symmetric, trace-free tensors to 1-forms. Since this latter operator 
is underdetermined elliptic, it is well-known that its adjoint is overdetermined elliptic and that the 
composition of these two operators as above is elliptic. 

These considerations lead to the following definition of the associated system, given here in 
index-free notation for ease of presentation: 

Ric"(S + h) ^ S(h,X,T) 

(17) 

divs+hoS{h,X,T)^0 

where S{h, X, T) will be called the York operator. As will be shown in due course, the map defined 

by 

$(/i, X, T) EE {mc"{5 + h)- S{h, X, T), div5+/i o S{h, X, T)) (18) 

on appropriate Banach spaces has a bounded, elliptic linearization in the h and X directions and 
as a result, the Implicit Function Theorem can be used to find solutions h{T) and X{T) as smooth 
functions of sufficiently small tensors T. 



3.3 Choosing the Banach Spaces 

Before proceeding with the solution of the equations (^^, it is necessary to specify in what Banach 
spaces of tensors the equations arc to be solved. The notion of asymptotic fiatness in R'^ should 
be encoded rigorously into the function spaces by requiring that the relevant objects belong to 
a space of tensors with built-in control at infinity. Furthermore, the spaces should be chosen to 
exploit the Fredholm properties of the operators appearing in the map $. Both these ends will be 
served by weighted Sobolev spaces, which are defined as follows. 

Let T be any tensor on R'^. (This tensor may be of any order — the norm || ■ || appearing in 
the following definition is then simply the norm on such tensors that is induced from the metric 
of R'"'.) The H^'P Sobolev norm of T is the quantity 

, 1/2 




where (j{x) = (1 -I- r^)^/^ is the weight function and = (x^)^ -I- {x^Y + {x'^Y' is the squared 
distance to the origin. Note that Bartnik's convention for describing the weighted spaces is being 
used (the reason for this is psychological: if / G H^^l^ and / is smooth enough to invoke the Sobolev 
Embedding Theorem (see below), then f{x) = o(r^) as r ^ oo, which is easy to remember — see 

for details). An appropriate choice of k and /3 for use in the Main Theorem will be made below. 

The space of i?'^-'' functions of R'^ will be denoted by H'^'^CR?) and the space of H'^'^ sections 
of a tensor bundle B over R^ will be denoted by H'''^{B). As an abbreviation, or where the 
context makes the bundle clear, such a space may be indicated simply by H^'^ . Note also that the 
following convention for integration will be used in the rest of this paper. An integral of the form 
/r3 /i in the definition above, denotes an integral of / with respect to the standard Euclidean 
volume form. Integrals of quantities with respect to the volume form of a different metric will be 
indicated explicitly, as, for example, J-^^ /dVolg. 

The spaces of H'^'^ tensors satisfy several important analytic properties and the reader is asked 
to consult Bartnik's paper, or others on the same topic [|[ |l^, 11 , for details. The three most 



important properties that will be used in the sequel are the Sobolev Embedding Theorem, the 
Poincare Inequality and Rellich's Lemma; these will be restated here for easy reference. 
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1. The Sobolev Embedding Theorem states that ii k > ^ and T is a tensor in H^'^ ^ then T i 



IS 



C". Furthermore, if the weighted C^' norm of a function / is given by 

fc 

ll/llc|-Ell^'/^"''^'llo' 

where ||/||o = sup{|/(x)| : a; G R^}, then in fact, T e and \\T\\co < C\\T\\h'^.p, 

2. The Poincare Inequahty states that if /3 < 0, then 

||/||//o,^ <C7||V/||^o.,-i, 
whenever / is a function in ^^''^(R^). 

3. The Rehich Lemma states that the inclusion H^'^{B) C (B), for any tensor bundle B, 
is compact when k' < k and (3' > (3. In other words, if Tj is a uniformly bounded sequence 
of tensors in H'^'^ , then there is a subsequence Tii converging to a tensor T in H'' . 

Remark: The constant C appearing in the estimates above is meant to depend only on the 
dimension n. In the remainder of this article, any constant depending only on n will be denoted 
by a generic C, unless it is important to emphasize otherwise. 

In addition to the three properties above, two important results that are valid in weighted 
Sobolev spaces will be needed in the sequel. The first concerns integration. 

Duality Lemma: If u G H'-'''{'R.'^) and v G ij'~2,-7-3^ then the integral u ■ v is well defined. 
Furthermore, the functional analytic dual space of H'^''^ (RP) is isomorphic to H'^' ''^'^(R'^) under 
the pairing v ^ 4>v where (j)v{u) — /j^3 u ■ v. 

Proof: Choose u and v as in the statement of the lemma. Then by Holder's inequality, 

\U-V\< [ |u|a-^-3/2 . |^|^-(-7-3)-3/2 

1/2 / /• \ 



R3 / \JR3 

< CXD . 

The product u ■ v is thus in and so its integral is well defined. The statement about duality 
follows from the Riesz Representation Theorem for and the inequality above. See ^ for 
details. □ 

The second result concerns the Fredholm properties of certain linear, elliptic partial differential 
operators on weighted Sobolev spaces. 

Invertibility Theorem: Suppose B is any tensor bundle over R'^ and let Q : H^'^{B) 
H^~^'^^^[B) he any linear, second order, elliptic, homogeneous, partial differential operator with 
constant coefficients mapping between weighted Sobolev spaces of sections of B, and k >2. Then 
Q is surjective i/ /? ^ Z and f3 > —1 and injective if (3 ^ 7i and (3 < 0. It is thus bijective when 
(3 G (—1,0). The operator Q is not Fredholm if (3 G 7i. 



Proof: The proof of this result can be found in |10 , but see also |^ for an excellent discussion of 
the intuitive foundation underlying the theory of elliptic operators on weighted spaces. □ 

Choice of Banach spaces 

Denote by iS'2(R'^) the symmetric tensors over R"^ and by Ai(R'^) the 1-forms of R"^. Solutions 
of the associated system will be found in the following Banach spaces. Pick any (3 G (—1,0) and 
any fc > 4; then 
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• metrics 6 + h will be found so that h G H''^'^ {S'^{R^)); 

• l-forms X will be found in H''-^'>^~^{A^CR^)); 

• tensors T will be found in H''-'^^'^-'^{S'^(R^)). 

The preceding choice of Banach spaces will be justified in the next section by showing that 
solutions of the associated system exist in these spaces. However, an argument can be made 
right now that suggests that the spaces above are indeed the correct ones in which to expect 
to find solutions. First, in order to ensure that the metric S + h is asymptotically flat, h must 
decay as r — > oo, and this holds by the Sobolev Embedding Theorem when f3 < 0. Next, a 
non-trivial, asymptotically flat metric satisfying the constraint equations must satisfy the Positive 
Mass Theorem and consequently must have non-zero ADM mass. Thus the term in the 
asymptotic expansion of h must be allowed to be non-zero, which by the Sobolev Embedding 
Theorem imposes the further requirement that (3 > —1. Furthermore, fc > 4 implies that the 
Sobolev Embedding Theorem can be applied to the second derivatives of the metric, and thus the 
curvature of the metric decays pointwise as r — » oo. Finally, the h, X and T quantities are chosen 
in different Sobolev spaces because of the differing numbers of derivatives taken on these quantities 
in the associated system. For instance, the reduced Ricci curvature operator is homogeneous and 
of degree two and thus sends a metric in H'^'^ to a tensor in ij'=-2,/3-2_ ^j^g operator S{h,X,T) 
is homogeneous but is only of degree one in X and of degree zero in T; it thus maps to i/'''-2,/3-2 
only when the weightings on X and T match together properly and match the weighting on the 
metric h as in the choice above. 

3.4 First Attempt to Solve the Associated System 

The Implicit Function Theorem, the tool which will be used to solve the associated system, is 
restated here for ease of reference. 

Implicit Function Theorem: Let ^ : A x B ^ C be a smooth map between Banach spaces 
and suppose that $(0,0) = 0. // the restricted linearized operator _D$(0, 0)|^^^q-j : A C is an 
isomorphism, then there exists an open set U G B containing and a smooth function (/) : U —i- A 
with 0(0) = so that $((/)(&), 6) = 0. 

For an excellent discussion and proof of this theorem, see Q . In order to use this theorem, let 

A = {{KX) e i7'"'^(S'2(R3)) X H''-^-'^'\A\R^))} 
B={TeH''-^^f'-^{S^{R'))} 

C = i/'^-2^'='-2(52(R3)) X ff'^-3.'3-3(Al(R3)) ; 

then the linearization of the operator <f> in the A direction at the origin must be calculated and its 
mapping properties understood. 

The linearization of $ is actually quite simple when evaluated at the origin because the covariant 
derivative of the Euclidean metric is trivial. The only nonlinearities in $ occur in the second order 
terms of the reduced Ricci operator and in terms that are quadratic in the derivatives of the 
metric (such as in products of Christoffel symbols or in the connection terms). It is thus easy 
to see that the linearization of a covariant derivative operator at the Euclidean metric is just the 
Euclidean derivative operator, and it is a straightforward matter to deduce from the definition of 
the associated system in ( p^ ) that the linearization of $ in the A x {0} direction is 

moAO)ih,x,o)^[-lfJll^^'^^), (19) 

where A is the Euclidean Laplacian and C is the Euclidean conformal Killing operator. 
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Denote by Ps the operator D$(0, 0, 0)(-, •, 0). It is a bounded linear operator between the 
appropriate weighted Sobolev spaces because of the way in which the weights were chosen in Section 
|3.3| . To determine whether Ps is an isomorphism, one appeals to the Invertibility Theorem. Recall 
that the weight /3 in the domain spaces of Ps has been chosen between —1 and 0. 

Injectivity of Ps 

Suppose {h,X) belong to the kernel of Ps{h,X). In other words, {h,X) solves the equation 
Psih,X) = {0,0),or 

-^Ah~jC{X)=0 
div o C{X) ^ . 

Since the operator divo£ : _ff'^~i''^^i(A^(R^)) i?'''^^'^^^(A^(R^)) is a linear, elliptic, homoge- 
neous, constant coefficient operator of second order, the Invertibility Theorem applies, and since 
/3 — 1 £ (—2, —1) when /? G (—1,0), it is thus injective. Hence X — 0. The remaining equation 
now reads Ah = and again, since A : ^'^■^(^^(R^)) ^ H''-^'^-^{S^{'R^)) and /3 G (-1,0), A 
is an isomorphism and thus h = 0. Hence Ps is injective. 

SURJECTIVITY OF Ps 

Although the operator Ps is injective, it is not surjective. First note that the Invertibility 
Theorem does not guarantee surjectivity in the same way that it guaranteed injectivity. To see 
this, attempt to solve the equations Ps{h,X) = {f,g) for any / G if'''~^'''~^(5^(R^)) and g G 
jjfe-3,/3-3|^^i^j^3-j^^ In other words, consider the system of equations 

-\Ah~C{X) = f 
div o C{X) = g . 

Because /3 — 1 G (—2,-1), the operator div o £ is not necessarily surjective according to the 
Invertibility Theorem. The full equations Ps{h,X) — {f,g) can thus not necessarily be solved. 

To show that Ps actually does fail to be surjective, it is necessary to show that the dimension 
of its cokernel in H'^''^ (S'^{'R^)) x H''^^''^^^ (A^(R^)) is strictly greater that zero. First, note that 
if Xg satisfies div o C[Xg) = g, then the remaining equation —^Ah = C{Xg) + / can be solved by 
the Invertibility Theorem since the weight /3 is chosen such that A is an isomorphism. Thus the 
dimension of the cokernel of Ps is equal to the dimension of the cokernel of div o £ as an operator 
between H^-^^I^~^{K^{H?)) and 7J*'-3'^-3(A1(R3)). 

To characterize the cokernel of div o C, one appeals to general, function-theoretic properties of 
linear, second order, homogeneous, elliptic operators on weighted Sobolev spaces. The following 
lemma and its proof show how this is done. 

Cokernel Lemma: Suppose B is any tensor bundle over R^ and let Q : H^''^{B) — *■ H^^^'^^^{B) 
be a linear, second order, homogeneous, elliptic operator mapping between weighted Sobolev spaces 
of sections of B where k >2 and 7 ^ Z, 7 < —1. The image of the operator Q is the space: 

Im(Q) = G iJ'^'"2''^"2(B) : J {w,z)^Q V z G Ker (Q*; -1 - 7))| , (20) 

where the inner product (•,•) is induced on B from the Euclidean metric o/R'^, the operator Q* 
is the formal adjoint of Q, and Ker(Q*; —1 — 7) is its kernel as an operator from H^'~^^^ {B) to 
jy'=-2.-3-T(B). 

Proof: Denote the space on the right hand side of equation (^o|) by W . Suppose that k = 2 and 
consider first the containment Im(Q) C W. Choose Q{y) G lm{Q) and z G Ker{Q*; —1 — 7). Since 
Q{y) G H^'^^^{B), the integral J-^^iQiy), z) is well defined by the Duality Lemma. 
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Claim: This integral equals {y , Q* (z)) . 

Proof: The equality of the integrals on smooth, compactly supported sections of B is true by 
definition of the adjoint. The equality of the integrals for H'^''^ sections follows because 
sections of B are dense in H''''* sections of B |5j. 

The integral {Q{y),z) is thus zero and so Q{y) G W. 

The reverse containment W C lm{Q) is proved as follows. Suppose wq belongs to W; thus, 
wq e H^'^^'^{B) and satisfies Jj^3{wo,z) — for all z G Ker{Q*; —1 — 7). Suppose also that 
Wo ^ Im((5). Since Q is elliptic, lm{Q) is closed; thus by the Hahn-Banach theorem, there exists a 
linear functional (j) on H^'''^^{B) so that 4>{wo) but 0|jjjj(q) = 0. Again by the Duality Lemma, 
there is a unique zq G H°'^^^'^{B) so that (f>{w) — /j^3(w, zq) for all w G H^'^~^{B). Therefore, 
'^lim(Q) ~ ^ implies that 

= <^(Q(y)) 

{za,Q{y)) 

(Q*(^o),2/) 



/R3 

for all y G H'^'^{B). Thus Q*{zo) = or G Ker(Q*;-l - 7). But now, the assumptions 
(t>{wo) ^ and J-^3{wq,z) = for all z G Ker{Q*; —1 — 7) are mutually contradictory. Thus it 
must be that wq G lm{Q). Finally, the extension to fc > 2 follows in a similar manner by standard 
functional analysis. □ 

Apply this theorem to the operator Q = div o C with 7 = /3 — 1. Now, Q* = Q, so in order to 
solve the equation div o C{X) — g, the tensors g must satisfy the constraints 



JR3 



0, 



where Y is any tensor in the kernel of the operator div o £ in the space H^^^'^^ ''(A^(R'^)). 

The kernel of divo£ is well known and consists of 1-forms dual to the the conformal Killing fields 
of R^. There are precisely ten linearly independent families of such vector fields: the translation 
vector fields, the rotation vector fields, the dilation field and three so-called special conformal 
Killing fields (these correspond to transformations of the form ioToi, where i is the inversion with 
respect to the unit circle and T is a translation) . The asymptotic behaviour of these vector fields 
can thus be computed exactly: the translations have constant norm, the rotations and dilations 
have norm growing linearly in the distance from the origin, and the special vector fields have 
quadratic growth in the distance from the origin. Since —1 — 7 G (0,1) when [3 G (—1,0), the 
only 1-forms dual to the conformal Killing fields in ff'^-^'-^ '(A^(R'^)) are thus those spanned 
by the translation 1-forms dx^, dx^ and dx^ . Consequently, the image of Q = div o C in the space 
^'£-3,7-2 ^yYi(R3)) can be characterized as follows: 

Im(divo/:) = jc, e i/'=-3,/3-3^^i(j^3)^ . J ^0, a= 1,2,3| , 

where ga are the components of g in the standard coordinates of R^. 

The conclusion that can be drawn from the analysis of this section is that the equation 
^{h,X,T) = (0,0) is not solvable near (0,0,0) using the Implicit Function Theorem. The non- 
surjectivity of the linearized operator at (0, 0, 0) is the essential obstruction. The best that can 
be achieved using the Implicit Function Theorem is thus that the equation ^{h, X, T) = (0, 0) can 
be solved up to a term that is transverse to the space /m(div o C). It will turn out that this is 
nevertheless sufficient for solving the full equations as a result of the compatibility conditions built 
into the equations. But in order to show this, the associated system defined in the previous section 
must be modified somewhat. 
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3.5 Reestablishing Surjectivity and Solving the Associated System 

In order to modify the associated system appropriately, first note that if'^~'^'''~'^(A^(R'^)) can be 
written as /m(div o C) (B W in many different ways; but in each case, is a three dimensional 
subspace of iJ'^~^'''^^(A^(R^)) whose members do not integrate to zero upon taking the Euclidean 
inner product with the translation 1-forms. One such choice is 

W = span {^dx'^l^^i 2,3 ' 

where <p is any smooth, positive function of compact support whose integral over R"^ is equal to 1. 

Again, denote the domain space of the operator $ by A. The previous paragraph suggests that 
one should attempt to construct a new associated operator that extends $ in such a way that 
$' : A X R^ — > Im{Ps) ® W, where the additional R'^ factor in the domain should map under the 
linearization D<i>' at the solution (0,0, 0;0) G A x R'^ onto the W factor in the image. If such a 
construction is possible, then the equation ^'{h, X,T; A) ~ (0,0) can be solved using the Implicit 
Function Theorem. 

Construct the operator : A x R'^ — > H''~^'^^'^ {^A^ (11^)) according to the prescription 

^\h,X,T;X) ^ (^Ric"{S + h) ~ S{h,X,T), divs+ho S{h,X,T) ^Y^^Xa(t>dx'''^ , (21) 

where, as before, Ric^ is the reduced Ricci operator and S'(-, •, •) is the York operator. The 
linearization of <&' at (0, 0, 0; 0) in the directions transverse to the T direction is easily seen to be 



D$'((5,0,0;0)(/i,X,0; A) = ^-^Ah- C{X), div o /:(X) - ^ A^^da;'^^ 



(22) 



Denote this new operator by P^. It is still bounded because has compact support, and it is now 
also bijective by the following arguments. 

Injectivity of Pg 

Suppose Pg{h, X; X) = (0,0). Integrate the components of the second equation; by the diver- 
gence theorem for the Euclidean metric (valid because constant functions can be integrated against 
jjk-3,p-3 functions when /? e (—1,0) according to the Duality Lemma), the divergence terms in- 
tegrate to zero, yielding Aa = for all a. The argument that both X and h are then equal to zero 



follows as in Section 3.4 



Surjectivity of Pg 



Suppose that Pg{h,X; A) ~ {f,g)- First choose the components Xa so that 

(.9a -I- Aa <?!>)= 



R3 

for each a. The equation div o C{X) = g — J2l=i ^a(t>dx'^ can then be solved for Xg according to 
the characterization of the image of the operator div o C from the previous section. The remaining 
equation —^Ah = 'C{Xg) + / can then be solved because /? G (—1,0) makes A an isomorphism. 

The Implicit Function Theorem can now be invoked to solve the equation $'(/i, X, T; A) = (0, 0) 
near (0, 0, 0; 0). To be precise, there is a neig hbourhoodZ/^ C H''-'^'^~'^{S'^{Ii^)) with the following 
property. If T G then there is a metric S + h{T) with h{T) G H''^'^ {S'^{R^)) , a covector field 
X{T) G i?'=~i'^~i(Ai(R^)), and three real numbers Xa{T) so that <i>' {h{T), X{T),T; X{T)) = 
(0,0). Furthermore, the various functions T i— > h{T), etc. are smooth in the appropriate Banach 
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space norms. In particular, there exists a constant C so that 



\\X\\uk-l.l3-l < C\\T\\ffk-2,l3-2 (23) 

||A||r3 <C||T||h.-2,^-2, 
where || • ||r3 denotes the standard Euchdcan norm of R"^, as long as T G iY. 

3.6 Satisfying the Harmonic Coordinate Condition 

Section |3.5| shows how the associated system (|l^) can be modified in such a way that it can be 
solved using the Implicit Function Theorem. This procedure results in a family of solutions of the 
equations 

Ric"{5 + h) = S(h,X,T) 

(24) 

div5+,i o S{h, X, T) = X(j), 



where A = X]a=i da;°. It remains to show whether the original equations (16) are satisfied by the 
solution S + h and S{h, X, T). This will be done by showing that the compatibility conditions built 
into the extended constraint equations (i. e. the Bianchi identity only, since the time-symmetric 
assumption has eliminated the other compatibility condition) actually ensure that if {h, X, T; A) 
solves (|2j), then A = and h + S satisfies the harmonic coordinate condition. Therefore Ric^{5 + 
h) = Ric{S + h) and solutions of ( |2^ ) are indeed solutions of the full equations. 

To prove this claim, assume instead that both A and the quantities are nonzero. Ar- 
gue towards a contradiction as follows. First, write g = S + h for short. The Bianchi identity 
divg (i?i!c((7) — ^R{g)g^ = 0, applied to equation (^) defining the reduced Ricci operator yields the 
identity 

— i^ab ~ 9ab) . = {j^ a;b + ^b-a ^ ^%hab) . 



which is equivalent to 



r,^,'^ + RtTa = 20Aa , (25) 



after using the modified associated system and commuting covariant derivatives appropriately. If 
Qh denotes the operator Ua i-^ /S.s+h'Ua + [Ric{5 + h)Wub, then ( p5| ) asserts that 2(j)Xa is in the 
image of H^-'^-^-'^{K^{n^)) under Qh, because h G H^^'^ {S'^{B?)) and the F" are obtained from 
5 + h hy differentiation. This, however, can be shown to violate the following basic result about 
elliptic operators. 

Stability Lemma: Let B he a tensor bundle over Y{? and let : H'''^{B) — * H'~^''*~^{B), 
e G [0,1], be a continuous family of linear, homogenous, second order, elliptic operators, for all 
7 < —1. Furthermore, suppose Qg is uniformly injective for any e whenever 7 < —1; i. e. for each 
7 ^ Z, 7 < ~1, there is a constant C independent of e so that \\Qi^(y)\\fji-2,-,-2 > C\\y\\Hk,-i ■ If 
z ^ Im((5o); then there exists Eq > so that z ^ YtcJ^Q^^ for all e < Sq. 

Proof: Suppose the contrary; then for some 7 < — 1, there exists a sequence — > and a sequence 
yi G H^^^{B) so that z = Qaiyi)- By the uniform injectivity of Q^, \\yi\\H'--' ^ C'|jz||^(-2,7-2 and 
is thus uniformly bounded. By Rellich's Lemma, there exists a subsequence yi' which converges 
to an element y in H^^^'^^p, where p is small enough so that 7 + p < —1. Again, by uniform 
injectivity. 
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< C\\{Q^' - Qy){yj:)\\Hi-2.-,+P-2 

< C\\Qt' — Qj'llop • llz/j' IIh'.^+p 

by the continuity of and the uniform boundedness of yi. Here, || • ||op denotes the relevant 
operator norm. The subsequence y^/ is thus Cauchy in the H^^^+p norm and so y.j' ^ y in this 
norm. But now, 

z = hm Qe,,{yi') = Qo{y) , 

i — *oo 

contradicting the fact that z Im(Qo)- □ 

In order to derive a contradiction from ( p5| ) using this lemma, the uniform injectivity of Qh 
must be established and it must be shown that (j)Xa does not belong to the image of Qo- 

Uniform Injectivity of Qh 

Suppose that gft(u) = Oforw e H''-'^^^ {A^{R^)) where 7 < -1. In other words, r^.^^+i^^r^ = 
0. From this, one easily deduces 

-Ag\\u\\^ = 2(i?,fcMV - \\Vuf) . (26) 

Before continuing, recall the following facts about Green's identity in weighted Sobolev spaces. 
If functions u and v are chosen such that v G H^'''(R.^) and u G ff'^'~^~^(R^) for some 7, then 
the integrals appearing Green's identity for a general metric g on a large ball Br, that is 

/ uAgt;dVolg+ / VM-VudVolg== / u|^dAg, (27) 

where dAg is the area form of the metric g, are all well defined as r — > 00. Thus by applying a 
density argument as in the proof of the Cokernel Lemma, one can conclude that 

wAgW dVolg + / Vm • Vw dVolg ^ , 

in the limit of ( p7| ) as r — > cxj. 

With this in mind, integrate both sides of equation (|2^) against the volume form of the metric 
g — 5 + h to obtain 

-\ [ Ag\\ufdY0\g= [ i^afcuVdVolg - / | | | | ^ d Volg . (28) 

^ JR3 JR3 JR3 

bmce u S 

Green's Identity can be applied to the left hand side of ( p8| ) when 1 G ^' ^. 
This is true since 7 < — 1; thus the integral of the left hand side of (pa) is zero. Consequently, 



0< / II i?Jc(.9) II IjufdVolg - / ||Vu||MVolg 

JR3 JR3 

</ ||ffic(.9)|| ||uf dVolg - C / II V||m|| f dVolg (29) 

JR3 JR3 

for some constant C, by the Cauchy- Schwartz inequality and straightforward algebra. Next, assume 
that h is small in a pointwise sense (this assumption follows from the Sobolev Embedding Theorem 
if h is sufficiently small in the H'^'l^ norm and k > |). In fact, assume that h is sufficiently close 
to so that all norms, derivatives and volume forms of the metric g can be replaced by their 
Euclidean counterparts (at the expense of changing C of course). Finally, since ||u|| is a scalar 
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function, the derivative operator in (|29| ) can be replaced by the Euchdean derivative operator 
without introducing lower order terms. Thus, there exists a new constant C so that the estimate 

0< / \\mc{g)\\\\ur^C I ||V||u|lf (30) 

holds, where the norms and derivatives appearing here are those of the Euclidean metric. Next, 
Ric((7) £ f{k-2,t3-2 i^gcause g — 6 & H'^'^ . But since fc > |, the Sobolev Embedding Theorem gives 
Ric(g) e That is, 

sup||Ric(g) •(j-'^+^ll < C < oo, 

which implies that 

sup ||Ric((7) ■ cr^ll < C < oo , 

R3 

since f3 < 0. Finally, apply the Poincare inequality for weighted Sobolev norms to the function ||u|| 
to deduce 



/ 

^R3 



\\mcig)\\\\ur<\\mcig)-a^o f Wufa 

JR3 



2„-2 



<q|Ric(g)||co / ||V||^.|| 

JR3 



<C\\g^S\\c2 I II V 

JR3 



M 



<c||;i||^.,. / ||v||z.|ir (31) 

JR3 



again by the Sobolev Embedding Theorem and the fact that /? < 0. Using (|3l|) in inequality 
leads to the contradiction because the preceding estimates imply 



Q<{C\\hU..,-l) / ||V||^.||||^ 

JR3 

while if ||ft,||^fc,f! is sufficiently small, the right hand side above is clearly negative. Avoiding 
this contradiction requires V||u|| = 0. But since the Sobolev Embedding Theorem applied to 
u G H^~^^^ shows that ||u|| decays at infinity when 7 < —1, it must be true that u — Q. 

The operator acting on H^^^''i 1-forms is injective for all 7 < —1 whenever h is sufficiently 
close to zero in the H^^l^ norm. The uniform injcctivity follows in the standard way from the 
injectivity of each Q/j and the fact that the constant in the elliptic estimate for these operators is 
independent of ft., again provided h is sufficiently near to 0. 

Image of Qq 

The (pX term in (^J) was specifically chosen in Section 3.5 to satisfy the integral condition 
J^3{\(j),dx^) ^ (since Aa = for all a). This condition ensures that indeed 2</)Aa is not in the 
image of the operator Qo = acting on the space of iJ''^ 1-forms of R'^ because the image of 
in H'"'^ for 7 < — 1 is perpendicular to the harmonic polynomials of degree less than the nearest 
integer less than 7, and this always includes the constants. 

The Stability Lemma thus applies to equation ( |2^ ) and implies that can not be in the image 
of Qh when h is sufficiently small in the H^^P norm, unless of course A = 0. Now, by the injectivity 
of the operator Qh, this in turn implies that ||r|| = 0, or that = for each a. Consequently, 
the harmonic coordinate condition for the metric 5 + h is satisfied, and as indicated earlier, this 
implies that the the metric 5 + h and the tensor S{h, X,T) satisfy the time-symmetric extended 
constraint equations. This completes the proof of the Main Theorem. □ 

Acknowledgements: I would like to thank Helmut Friedrich for suggesting the topic of the 
conformal constraint equations to me as well as for providing me with guidance and insight. I 
would also like to thank Piotr Chrusciel, Jorg Frauendiener and Niall O'Murchadha for their 
valuable comments and suggestions. 



20 



References 



R. Abraham, J. E. Marsden, and T. Ratiu, Manifolds, Tensor Analysis, and Applications, 
second cd., Springer- Vcrlag, New York, 1988. 

L. Andersson and P. Chrusciel, On "hyperboloidal" Cauchy data for vacuum Einstein equations 
and obstructions to smoothness of Scri, Commun. Math. Phys. 355 (1994), 1-102. 

, Solutions of the constraint equations in general relativity satisfying "hyperboloidal 



boundary conditions", Dissertationes Mathematicae 161 (1996), 533-568. 

L. Andersson, P. Chrusciel, and H. Friedrich, On the regularity of solutions to the Yamabe 
equation and the existence of smooth hyperboloidal initial data for Einstein's field equations, 
Commun. Math. Phys. 149 (1992), 587-612. 

Robert Bartnik, The mass of an asymptotically flat manifold, Comm. Pure Appl. Math. 39 
(1986), 661-693. 

William M. Boothby, An Introduction to Differentiable Manifolds and Riemannian Geometry, 
2nd ed.. Academic Press, Inc, 1986. 

Adrian Butscher, Perturbative asymptotically fiat solutions of the extended constraint equa- 
tions in General Relativity, In preparation. 

Murray Cantor, Elliptic operators and the decomposition of tensor fields, Bull. Amer. Math. 
Soc. 5 (1981), no. 3, 235-262. 

Jian Guo Cao and Dennis DeTurck, Prescribing Ricci curvature on open surfaces, Hokkaido 
Math. J. 20 (1991), no. 2, 265-278. 

Y. Choquet-Brouhat and D. Christodoulou, Elliptic systems in Hg^g spaces on manifolds which 
are Euclidean at infinity. Acta Math. 146 (1981), 129-150. 

D. Christodoulou and N. O'Murchadha, The boost problem in general relativity, Commun. 
Math. Phys. 80 (1981), 271-300. 

Dennis DeTurck and Hubert Goldschmidt, Aspects of prescribing Ricci curvature, CR- 
Geometry and Overdetermined Systems (Osaka, 1994), Math. Soc. Japan, Tokyo, 1997, 
pp. 60-84. 

, Metrics with prescribed Ricci curvature of constant rank. I. The integrable case, Adv. 



Math. 145 (1999), no. 1, 1-97. 

Dennis M. DeTurck, Metrics with prescribed Ricci curvature. Seminar on Differential Geome- 
try, Princeton Univ. Press, Princeton, N.J., 1982, pp. 525-537. 

Jorg Prauendiener, Gonformal infinity. Living Rev. Relativity 3 (2000), no. 4, Online 
Article: cited in December 2001, http://www.livingreviews.org/Articles/Volume3/2000- 

4frauendiener / . 

Helmut Friedrich, Cauchy problems for the conformal vacuum field equations in general rela- 
tivity, Commun. Math. Phys. 91 (1983), 445-472. 

, Asymptotic structure of space time. Recent Advances in General Relativity (A. Janis 

and J. Porter, eds.), Einstein Studies, vol. 4, Birkauser, 1992. 

, Einstein's equations and conformal structure, The Geometric Universe. Science, Ge- 
ometry and the work of Roger Penrose (S. Huggett, L. Mason, K. P. Tod, S. Tsou, and N. M. J. 
Woodhouse, eds.), Oxford University Press, 1998. 



21 



[19] , On the existence of n-geodesically complete or future complete solutions of Einstein's 

field equations with smooth asymptotic structure, J. Geom. Phys. 24 (1998), 83-163. 

[20] Emmanuel Hebey, Sobolev Spaces on Riemannian Manifolds, Springer- Verlag, Berlin, 1996. 

[21] Robert C. MeOwen, The behavior of the Laplacian on weighted Sobolev spaces. Comm. Pure 
Appl. Math. 32 (1979), 783-795. 

[22] Roger Penrose, Asymptotic properties of fields and spacetimes, Phys. Rev. Lett. 66 (1963), 
no. 10. 

[23] , Zero rest-mass fields including gravitation: asymptotic behaviour, Proc. R. Soc. Lond. 

A284 (1965), 159-203. 

[24] R. Schoen and S.-T. Yau, On the proof of the positive mass conjecture in general relativity, 
Commun. Math. Phys. 65 (1979), 45-76. 

[25] Giinter Schwarz, Hodge Decomposition — A Method for Solving Boundary Value Problems, 
Springer- Verlag, Berlin, 1995. 

[26] J. York, Gravitational degrees of freedom and the initial-value problem, Phys. Rev. Lett. 26 
(1971), 1656-1658. 



22 



